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We propose massive gravity as a holographic framework for describing a class of strongly interact- 
ing quantum field theories with broken translational symmetry. Bulk gravitons are assumed to have 
a Lorentz-breaking mass term as a substitute for spatial inhomogeneities. This breaks momentum- 
conservation in the boundary field theory. At finite chemical potential, the gravity duals are charged 
black holes in asymptotically anti-de Sitter spacetime. The conductivity in these systems generally 
exhibits a Drude peak that approaches a delta function in the massless gravity limit. Furthermore, 
the optical conductivity shows an emergent scaling law: jo-(cj)| ~ + B . This result is consistent 
with that found earlier by Horowitz, Santos, and Tong who introduced an explicit inhomogeneous 
lattice into the system. 



I. INTRODUCTION 

In systems with perfect translational symmetry, the 
particles cannot dissipate their momentum. Conse- 
quently, in the presence of a finite density of charge car- 
riers, there is a delta function in the AC conductivity 
at zero frequency. The DC conductivity is therefore in- 
finitely large. This unwanted result can be avoided by 
treating the charge carriers in the probe limit (i.e. as 
a small part in a larger system of neutral fields where 
they can dump momentum), or by introducing spatial in- 
homogeneities thereby breaking translational invariance 
explicitly. 

Once momentum dissipation has been introduced into 
the system, the results will be finite. This can be demon- 
strated by the Drude model of conductivity: a phe- 
nomenological theory that treats the charge carriers as 
classical particles which can bounce off a background ion 
lattice. The equation of motion. 
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where r is the mean free time between collisions, q is 
the electron's charge, E is the background electric field 
driving the current. The DC conductivity is then finite 



The transport properties of a large class of metals are 
reasonably well described by the Drude model. There 
are materials, however, whose optical conductivity de- 
viates from the simple Drude formula. In one class of 
high temperature superconductors, for instance, the ob- 
served conductivity in the normal phase follows a power 
law |cr(w)| cx (— zcj)"^/^ (see [I1I2])- These systems are 
strongly coupled and there is no simple quasiparticle- 
based Fermi liquid description. 

In the past fifteen years, there has been much progress 
in understanding certain strongly interacting quantum 
field theories in the 't Hooft limit using the AdS/CFT 
correspondence [3H5]. Translational symmetry breaking 
has been studied and Drude(-like) peaks were discovered 
[51-fTU]. Recent numerical calculations in holographic lat- 
tice systems gave evidence for universal non-Drude fre- 
quency scalings [TTl [T^] . 

Holography in itself is not doing any coarse-graining, 
therefore such calculations on inhomogeneous back- 
grounds require the solution of partial differential equa- 
tions^. This motivates the main goal of this paper which 
is to build a framework for translational symmetry break- 
ing and momentum dissipation in holography without the 
need for complex numerical computations. 
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II. HOLOGRAPHIC MATTER 



where n is the number density and m* the effective elec- 
tron mass. In order to compute the AC conductivity, 
one turns on a time-dependent electric field with angular 
frequency uj. This yields 
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The pole is shifted to the lower half plane and the DC 
conductivity is finite. 



In order to describe the holographic dual of strongly 
coupled matter, we are going to use a set of minimal 
ingredients, namely, Einstein-Hilbert action with a gauge 
field and a (negative) cosmological constant 
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^ For conductivity calculations on a homogeneous (Bianchi VII) 
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The equations of motion are solved by the following AdS- 
Reissner-Nordstrom geometry [18] 



dr^ 



^f{r)dt^ + dx^ + dy^ 



A{r) ^ —Jdt 

where the emblackening factor is 
/(r) = 1 - Mr^ + 
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This is a charged black brane with a horizon at rh which 
is the smallest positive root of f{r). Using the AdS/CFT 
dictionary, fi is identified with the chemical potential of 
the boundary theory. In the zero-temperature limit, the 
near-horizon metric becomes AdS2 x which governs 
much of the low-energy physics. 

This geometry describes a translationally invariant 
state. The electrical conductivity in the boundary theory 
can be computed by means of the Kubo formula 
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where (. . .) denotes the retarded Green's function. The 
currents are dual to the bulk gauge field. Thus, we need 
to compute the boundary-to-boundary two-point func- 
tion of the gauge field ax at finite frequency and zero 
spatial momentum. At non-zero charge density, due to 
the background gauge field profile, ax{r) mixes with the 
graviton gtx{^)- Let us rescale the variables such that 
rh = To first order in the perturbations, the Maxwell 
equation is 

2uj'ax + f [-KiJtY + 2 {a'xf + fO] - (9) 



where gf = gtx9^^ , and at the linear level g^ 
r — X component of Einstein's equations 



A^rax 
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The t — x component of Einstein's equations follows from 
the r — X equation above (by taking the derivative). 
We impose ingoing boundary conditions at the horizon. 
This corresponds to computing retarded Green's func- 
tions [19.. The two coupled differential equations can 
then be solved (e.g. by the method of matched asymp- 
totic expansions or numerically). In the ultraviolet re- 
gion, ax ^ a+ + a-T + 0{r^), and the Green's function 
is given by the ratio G = ^ . 

The result contains a Dirac-delta function in the real 
part of the conductivity. This is a consequence of the 
Ward identities for the translational symmetry. In the 
following, we resolve this delta function by letting mo- 
mentum dissipate. 

For more information on holographic matter, see the 
review papers 



III. BREAKING TRANSLATIONAL 
INVARIANCE 

In the charged black brane background from the previ- 
ous section, Ward identities for translational invariance 
in the x direction imply a shift symmetry in the gf field. 
This is why only derivatives of gf arose in the equations 
of motion. In order to to break translational symmetry, 
the shift symmetry must be broken. The simplest option 
is to add a mass term for the graviton. 
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where indices are raised using the diagonal background 
metric. Since the background is diagonal, Sgtx = gtx- 
The graviton mass term produces a linear gtx term in 
the t — x component of Einstein's equations. However, 
the r — x component is unchanged and the two equations 
are now incompatible unless the grx graviton component 
is also non-vanishing. This component carries an extra 
degree of freedom.^ 



A. Non-linear massive gravity 

If we add generic mass terms for the gravitons on a 
given background, then the theory will be plagued by var- 
ious instabilities, sometimes at the non-linear level. Re- 
cently, the authors of [26] constructed a theory where the 
Boulware-Deser ghost [57] was eliminated by introducing 
higher order interaction terms into the Lagrangian. (See 
also earlier works [55H52] . For a recent review of massive 
gravity in this context, see |33j.) 

In 3-1-1 dimensions, the theory has two dimensionless 
parameters, and it also depends on a fixed rank- 2 sym- 
metric tensor /, the reference metric. The usual dynam- 
ical metric will be denoted by g. For our purposes, we 
include a cosmological constant and a Maxwell field, 
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where Ci are constants, Ui are symmetric polynomials of 



^ An interpretation of this field is the following. In a semi- 
holographic effective theory I23H25| . one separates the bulk 



spacetime into a UV and an IR region, 

S = SuY{Gij,gjj) + SiB.{gij) 



(12) 



where G/j is the UV boundary value of the metric, and gjj is 
its value at a fixed intermediate cutoff scale, and I, J a {t, x, y}. 
The action is invariant under coordinate transformations that 
change either g or G. In the low energy theory the two symmetry 
groups are broken down to the diagonal. Finally, the Goldstone 
bosons corresponding to the broken axial symmetry are the radial 
integrals of the g''^ fields in the UV region. 
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the eigenvalues of the 4x4 matrix /C^i^ = Vs^^/qi^ 

[K] (14) 

Ui ^ [/C]4 - &[JC^][JCf + 8[/C3][/C] + 3[/C2]2 _ 6[/C^] 

The square root in K, is understood to denote matrix 
square root, i.e. {^/AY ^^{'s/AY ^ = A^^. Rectangular 
brackets denote traces: [K] = K,^ ^. As to —j' 0, we re- 
cover the translational invariant action in section Ull 

If the reference metric is flat, we can express it via a 
coordinate transformation (f>°'{x) using r/ab, 



(15) 



Different choices for the 0" Stiickelberg fields correspond 
to different gauges. The unitary (or physical) gauge is 
defined simply by 0° = x'^S'^. 

In this paper we will be interested in the case of a 
spatial reference metric (in the basis {t,r,x,y)) 



= ifsp)fi^ = diag(0,0, 1, 1) 
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Note that the action remains finite since it only contains 
non- negative powers of /^i/. The reason for using this 
singular metric becomes clear if we perform a coordinate 
transformation <f>°'{x) which preserves the x and y coor- 
dinates. 



6*-'- = 0*''-(t,r) 



(17) 



The reference metric and the action stay the same. This 
means that the spatial graviton mass term m?U{g, fsp) 
preserves general covariance in the t and r coordinates, 
but breaks it in the two spatial dimensions. This is ex- 
actly what we need. 

Since the reference metric is spatial, we only need two 
Stiickelberg fields cf)^ and 0^. They can be thought of 
as maps used in the Lagrangian representation of the 
degrees of freedom in a solid. In analogy to crystals, these 
dofs may be called ions. Perturbations of the Stiickelberg 
fields are the phonons 
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In this interpretation, the bulk is filled with a homoge- 
neous solid [S^ . Due to a gauge symmetry we can either 
set gri = and have tt* 7^ or vice-versa. 

In the ADM formulation [35], the metric is 
parametrized in the following way: N = (—5°°)"^^^, 
= .9oii 7ii — dij- Furthermore, let us define j'^^'jjk = 
61 and iV' = ^^^Nj. In terms of these variables, the 
spatial graviton mass term assumes the explicit form'^ 



sp. 
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V2 = Vdet(7-V)Vl-^^J- 



where 7 ^ = 'J^^^) , ^ = ^) \ and the spatial 
submatrix of the reference metric was kept in a general 
form: / = (t'"' l^") . a and P are two parameters (equal 

to ci and c^ in eqn. ( |13| ). 

Both V1.2 are invariant under spatial rotations. Note 
that they do not contain and thus the corresponding 
constraint in the ADM formulation is preserved. Note 
that the g'"* components do not show up in these terms 
either'* . 

Instead of using /C, we could consider mass terms made 
out of YZp" = a/ g/ia/"" where the inverse reference met- 
ric is set to f'^'^ = diag(0, 0, 1, 1). It is easy to check that 
the new functions Ui{JC) do not contain the lapse and 
shift fields. However, unlike Ui{JC), they are functions of 
■ We will not consider this option here. 

In this paper, we will not consider the delicate ques- 
tion of ghosts and tachyons. These investigations typi- 
cally depend on the background metric and the choice of 
parameters (see e.g. [3^]). We just note here that the 
Hamiltonian constraint can be preserved by redefining 
the fields — >■ n* using a transformation that is linear 
in N 



(20) 



and then choosing such that y/^U oc NU{n^,N,'^) 
is linear in N. For instance, assuming / — Id2x2, for 
V2 this can be done by setting = and — 
+ {n^Y^ det(7)~-'-/^-y7ijn*ri^ A^. Then, the mass term 
becomes linear in N 

NV2 = det(7)-i/y7,,n*nJ- -t- (^7"^ " ^""^) ^ (21) 

Thus, the Boulware-Deser ghost is eliminated from the 
theory. For more on this, see 



IV. GRAVITY BACKGROUND 



In the following, we study the massive gravity action 



(131 with the explicit mass term (19) and set 



= diag(0,0,F2,F2 



(22) 



(19) We will be looking for charged black brane solutions. 



^ There are only two terms, since the spatial gauge only allows for 
two non-zero eigenvalues for the matrix K. 



. . . even though A^' do appear. This is due to the asymmetric 
parametrization of the metric w.r.t. the r and t coordinates. 
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We obtain the following solutions to the equations of 
motion^ 



_drf_ 



Mr) = A* 1 - 



-f{r)dt^ + dx^ + dy' 



dt 



(23) 



(24) 



where the emblackening factor is 



fir) 



aF r + PF^m^r^ - Mr^ + A^r"* (25) 

2 4r^ 



For the equations of motion, see Appendix A. The hori- 
zon is located at where both /(r) and A{r) vanish. 

There are two dimensionless parameters: a and (3. F 
and m are redundant parameters and are only included 
for convenience. Note that if m = or F = 0, then the 
solution reduces to the AdS'-Reissner-Nordstrom solution 
of massless Einstein-Maxwell theory. The temperature is 
given by 



T = 



1 



3- 



(^^y + Fr^m^aL + (^Fr^)\ (26) 



Note that whenever /i < 2mFy/]3, the function T{rh) has 
a minimum at 



V3 



m^F^l3^l^ 



(27) 



and there is a corresponding minimal temperature T = 
r(rmin)- Below the critical size the black brane is unsta- 
ble. 

The geometry describes a finite density state with the 
entropy, charge and energy densities respectively given 
by 



4L^ 
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( 1 , Fm^rh(La+2FrHl3) 
1 + io^^2 



) 



i^Mrl + FLm'a 

h \ 



4F^m'^rh{La+2FrH/3y^ 
2Mrl-Fm'^rh{La+4Frhl3)+4: 



These quantities were obtained by computing the action 
for the (Euclidean) background with the UV divergences 



^ Note that using the reference metric, a new invariant can be de- 
fined: I"''' = gi^''d^(t>°d,,(p''. In unitary gauge, I^'' = gi^^S^S^, 
which is singular if the inverse metric is divergent. This pre- 
sumably leads to perturbative instabilities. Hence, in order to 
describe a black hole, one typically needs an ansatz for the ge- 
ometry in which the metric has no horizon singularities. On the 
other hand, in our spatial gauge does not depend on g*' since 
ft/j, = 0. Thus, wc will be able to use simple coordinate systems. 
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FIG. 1: Stability in parameter space. We set ~ L — m = 
F = 1 for the plot. Above the "wall of stability" P = - 
the entropy density is larger than the usual value ( "S = A/4" ] 
and numerical results indicate an instability. On the line /? = 

3 + FLm^rf^c 



2 '2 



, the maximal value of the chemical potential is 

zero. Between these two lines (yellow region) the system can 
be stable. The lines cross at (a, /3) — (~ „. ^, — , ^, , ). 

Beyond this point there may still be stable points. 



removed. This defines the grand canonical ensemble from 
which we get the above results (see [IS] for similar cal- 
culations on Ads'- Reissner- Nordstrom backgrounds) . By 
construction, e, s and p satisfy the first law of thermody- 
namics 



de = Tds + ^dp 



(28) 



Interestingly, the entropy density differs from its usual 



value So 
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unless 
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In this case, the energy and charge densities also simplify. 
This line will be called the wall of stability for reasons 
that will become clear later. 

Let us fix r/j = 1. At fixed graviton masses, p is max- 
imized if we set T = 0. On the line 



3 + FLm^rhOL 
F^m^rf^ 



(30) 



the maximal value of the chemical potential is zero (here 
p = as well). This will be called the "/x — 0" line. 
The ground state entropy is found to be 



s{T = 0) 



iL^n / m2 {FLr^a + 2F^rll3) ' 
Tf^ I ^ 2L{F'm?rha + &/L) 



(31) 



On the = 0" line it is equal to so/2. It would be 
interesting to find an interpretation of these results. 
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V. CONDUCTIVITY 

In order to compute the conductivity, we perturb the 
background^ 

ds^ ^ ds^ + 5t,(r)e^"* + g„(r)e*"* (34) 

A{r) ^ A{r) + a^{r)e"^'dx (35) 

The equations of motion are presented in Appendix B. 
The graviton mass does not appear in the Maxwell equa- 
tion. Its effects are communicated to the gauge field only 
through the coupling to the graviton fields. 

From the equations we can read off the masses of the 
grx and gtx fields. They are equal and depend on the 
radial direction 

rF 

rn^{r) — — - {aL + f3Fr) w? 



2L2 



Using this formula the temperature can be rewritten as 

At T = 0, the scaling dimension of in the infrared 
^^5*2 is given by 



^)' + 2LW(r,)) 



(36) 



A = 




16 (6 + Lm?rha) 
\2 + m?rh{?,La + 2rhli) 



17- 



(37) 



From this we get A = 2 (the rn = result) only on the 
wall of stability where /? = — . On the "/i = 0" line 
the formula gives A = 1. 

We can eliminate gtx from the equations and obtain 
two coupled second order equations for g^x and ax ■ These 
two equations have been used for the numerical calcula- 
tions in this paper. At the horizon, we impose infalling 
boundary conditions 



ax{r), grx{r) oc (r,,. - r) ^-"r 



(38) 



We set normalizable UV boundary conditions for the grx 
field^. This determines the wavefunctions up to a con- 
stant factor. We proceed to read off the Green's func- 
tion: ax ^ a+ + a_r -|- 0{r^) near the boundary, and 
then G{uj) = ^ as earlier. Finally, the Kubo formula 
gives the conductivity: = G{oj)/{iu}). 

In the general case, the conductivity exhibits a Drude 
peak as seen in FIG. [2j The size of the peak grows as 
m decreases. In the m — >■ limit, we recover the delta 
function at u — Q. 



® Equivalently, one can also consider the perturbation 



A(r) ~¥ A(r) + a:,{r)e^^'^dx 



(32) 
(33) 



b 

a 
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FIG. 2: Drude peak in the conductivity. The real and imag- 
inary parts are drawn in blue and purple, respectively. At 
larger frequencies, the conductivity approaches a constant. 
The parameters were set toa = — 1, /3 = 0, /i = 1.724, 
T = 0.1, m = 1, L = 1. 



A. Stability 

In this paper we will not attempt to thoroughly study 
the conditions of stability. There are certainly inconsis- 
tent regions in the parameter space, where the retarded 
gauge field correlator contains poles on the upper half- 
plane of complex frequencies. Numerical results indicate 
that this happens above the wall of stability (see FIG. 
IT]). Between the wall of stability and the "/j, = 0" line 
(the yellow region in FIG. [T]) the system may be stable. 



B. Emergent non-Drude scaling 



The optical conductivity differs from the simple 
|CT(a;)| oc u)~^ that is predicted by the Drude theory. In 
an intermediate regime T < ui < fi, we see an approxi- 
mate behavior best described by 



A 

uj'y 



B 



(39) 



where 7, A and B are 0(1) constants that depend on the 
a and /3 parameters. Numerical results gave _B < in all 
cases. See FIGs. 3(a)[ |3(b")| |3(c) for a sample numerical 
solution. In these figures, we have tuned the graviton 
masses such that 7 ~ —2/3. The power law behavior 
extends to larger and larger regions as the temperature 
is lowered. By changing the graviton masses, we obtain 
power laws with different exponents. As m — >■ 0, the peak 
becomes more and more Drude-like (i.e. 7 = 1). These 
results are very similar to those in [TI] [T^] . 



If we intend to use phonon fields instead of the grx graviton, 
then we may set an (equivalent) Dirichlet boundary condition 
on 7r'(r) at the UV boundary. (Due to a shift symmetry, 7r(r) 
itself does not appear in the equations of motion.) 
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FIG. 3: Non-Drude optical conductivity. There is an approximate power-law: ~ -A^ + B. The mass is tuned to 

L^m?a — —0.75 (and /3 = 0) so that the exponent 7 ~ 2/3 with an ofTset B ~ —1.2. The constants 7, A and B depend on the 
two param eters a and /3. Fig. 3(a) The blue and purple lines are the real and imaginary parts o f the conductivity, respectively. 

Phase of cr(a;). If B were 



Fig. 



3(b) 



The plot shows 



which gives the exponent if there is indeed a power law. Fig. 



. d(logi^) 

zero, then it would exactly be 60° due to causality and time-reversal invariance [2]. 



3(c) 



VI. DISCUSSION 

In this paper, we studied massive gravity as a holo- 
graphic framework for translational symmetry breaking 
and momentum dissipation. Instead of directly dealing 
with inhomogeneous fields (e.g. the metric) and their 
perturbations, we considered averaged quantities that 
satisfy 'renormalized' equations of motion. Ideally, fig- 
uring out what these modified equations are would be 
done by integrating out high-wavelength modes in the 
theory. However, that is a hard problem and instead we 
made a step by guessing their form by considering the 
symmetries of the system. 

We arrived at a holographic theory of solids based 
on Lorentz-breaking graviton mass terms. We com- 
puted conductivities in different cases. The conductivity 
showed a Drude peak at zero frequency. We also observed 
non-Drude power-laws in the absolute value of the opti- 
cal conductivity. These fat tails extended to frequencies 
comparable to the chemical potential. This ultraviolet 
effect seems to be unrelated to the physics that governs 
the DC conductivity, which in our model ultimately de- 
pended on the graviton mass only (at relatively small 
temperatures). It would be interesting to consider more 
general models which allow for a temperature-dependent 
DC conductivity. In order to do this, the reference met- 
ric may be promoted to a dynamical quantity (see related 
work f4(T ) . For instance, if we consider the simple ansatz 

-diag (0,0, i^(r),F(r)) (40) 

then the equation of motion for F{r) gives 

F(r)(X5..(r) = ^ (41) 

This corresponds to a constant m^{r) which is further 
equivalent to a shift in the cosmological constant. One 
might also add a kinetic term and a mass term for F{r) so 



that the r-dependence changes and it emulates the finite- 
momentum 'master field' in [in] . It would be very impor- 
tant to develop a quantitative correspondence between 
lattice perturbations and graviton masses (and perhaps 
higher order corrections in the action). 

We emphasize that we did not attempt to character- 
ize the instabilities in these systems. There are certainly 
inconsistent regions in the parameter space, where the re- 
tarded gauge field correlator showed poles on the upper 
half-plane. It would be important to understand under 
what conditions can the ghosts and tachyons be elim- 
inated from the theory. Such investigations generally 
depend on the background. Some instabilities are pre- 
sumably associated with the growth of inhoniogeneities 
with time (structure formation). 

Instabilities generically lead to other phases. A very 
interesting application of the results would be to study 
charged condensates as a model for supersolids or per- 
haps a pseudogap phase. It would also be interesting 
to find striped, dielectric, or insulating phases by chang- 
ing the reference metric, or study electron stars |41) in 
this context. One might also wonder whether there is a 
holographic analog of Cooper pairing and study how the 
effective phonon coupling changes in the radial direction. 

One can try to extend the theory to other dimensions. 
In three spatial dimensions one might use homogeneous 
Bianchi metrics as a spatial reference metric which ex- 
tends the number of possibilities. (For Bianchi spaces 
and holography, see e.g. [TOIH^HIS] .) 

In the paper, general covariance was only broken in 
the spatial dimensions. It would be extremely valuable 
to develop a theory where time translations are also bro- 
ken. This may be a first step toward constructing a holo- 
graphic model of Kolmogorov's 1941 theory of fully de- 
veloped turbulence. 

Computations in holography might look complicated 
at first. The physics is often elucidated by a semi- 
holographic approach |23| . It would be very useful to 
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develop an effective tfieory along the lines of [23]. At 
the technical level, analytical results would be extremely 
useful, perhaps by using matched asymptotic expansions. 
Finally, it would also be interesting to see if there are any 
relations to other, non-relativistic, holographic theories 
[46] . In particular, how the conductivity calculations on 
Bianchi VII spaces that produce a Drude peak |13.^ can be 
reformulated in our framework. This would presumably 
shed more light on the criteria of consistency in massive 
gravity theories. 

As we have seen in section |IV[ thermodynamical state 
variables are modified by the finite graviton mass. In 
particular, entropy does not generically follow the usual 
S — A/A law. It would be interesting to interpret these 
results and compute other related quantities (e.g. entan- 
glement entropies). 

In this paper, we substituted graviton mass terms for 
spatial inhomogeneities in asymptotically anti-de Sitter 
spacetime. As mentioned earlier in this section, gravi- 
ton masses may cause a shift in the value of the effective 
cosmological constant that is seen by perturbations. The 
relation of momentum dissipation and an effective cos- 
mological constant can be demonstrated in a more direct 
way as follows. In the absence of external forces, we can 
write the algebra corresponding to the Drude model as 



dt T 



{P^,PJ}=0 



where H is the Hamiltonian and {•, •} is the Poisson 
bracket. In 2-1-1 dimensions, the algebra spanned by 
{H,Px,Py} corresponds to Bianchi type V spaces where 



time plays the role of one of the three Bianchi dimensions. 
The simplest example for a spacetime whose Killing vec- 
tors obey this algebra is de Sitter space 



-"^dx^ 



(42) 



and a corresponding cosmological constant is given by 
A = 3/(2t^). We see two dual pictures emerging. We 
either have flat space with momentum dissipation, or de 
Sitter with conserved momentum. In the latter, we have 
traded momentum dissipation for the expansion of space: 
momentum is simply being inflated away. 

Using these ideas, one can calculate a 'mean free path' 
corresponding to the cosmological constant of our Uni- 
verse. We get A = cr w 3.4 Gpc. Amusingly, this is 
only a few dozen times bigger than the 100 Mpc 'lattice 
spacing' of the large-scale structures (above which the 
Universe is approximately homogeneous and isotropic). 
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Appendix A: Equations for the background 



Einstein's equations are supplemented by the graviton mass term, 



R 



Rt.. - + Ag^. -I- F^^F^ + ^F^f.F'^^ + m^X^, = 



where A ^ and 



X., 



(Al) 



(A2) 



with IC^^{g, /) = {.\/ f)^ Indices are lowered and raised by the metric g. The t — t component gives a differential 
equation for the emblackening factor 



2 4 

2r/'(r) - 6/(r) + 2aFLm^r + 2l3F^m^r^ - +6 = 



whose solution for /(r) is the one given in Section IV 



(A3) 



Appendix B: Equations for perturbations 



At the linear level, we obtain the following three equations 

2L\hL0^a.Ar) + fir) {2L\U{r)al{r) + 2L\ha'^{r)f {r) + itir^ugr,{r) - iir'g'.Ar) - 2ixvgt^{r)) = 



8 



gtx{r) {-2r^rlf"{r) + 8rrlf{r) - 16rlf{r) + AaFLm'^rrl + 4:/3F^m.^r^rl + fi^r^ + Url) + 
+2rrhf{r) {-2iiL\a'^{r) + m {-irujg'^,{r) + rg'l^{r) + 2g'^^{r))) = 

9rx{r) (/(r) {-2rl {r'^f"{r) - 4r/'(r) + 6/(r)) + 4rl {aFLm.^r + pF^m^r'^ + 3) + /x^r^) + 2r^rloj^) + 

+2irrlu) {rgl^{r) + 2gt^{r)) - AtnL^r^rhUJa^ir) = 

Note that the equations become real if we multiply grx by i. If m = 0, then we can consistently set g^x = and then 
the second equation becomes dependent on the third one. Since m > only introduces gtx, but not its derivatives, 
we can express gtx using the other variables. 
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